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Abstract. In this paper, we show that for an n- vertex graph G of genus 
g, the edge expansion of G can be determined in time n°' 9 '. We show 
that the same is true for various other similar measures of edge connec- 
tivity. 

1 Introduction 

1.1 Background and Motivation 

Edge expansion (known also as the minimum cut quotient, the isoperimetric 
number, or the flux of a graph) is a well-studied notion in graph theory and arises 
in several contexts of discrete mathematics and theoretical computer science. 
These include the explicit construction of expander graphs, the analysis of certain 
randomised algorithms, and graph partitioning problems. In this paper, we are 
concerned with giving an exact algorithm for determining the edge expansion 
(and other similar measures) of graphs embedded on surfaces. 

Throughout, we use the term graph to mean multigraph without loops, unless 
otherwise stated. For a graph G = (V, E) and e g E, we write e = ah to mean 
that the vertices a and b are the end points of e. 

For S a nonempty proper subset of V and S its complement, we define 

[S, S] G = {eeE: e = ab, a e S, b e S}, 

which we call an edge- cut of G. (The subscript is dropped when it clear which 
graph we are referring to.) For a cut [S,S] of a graph G = (V,E), define the 
balance of the cut to be b(S, S) :— min(|5|, |5|)/|V|. Note that the balance of 
a cut is a real number in the interval (0, h]. Two well-known graph cut prob- 
lems, which take into account the balance of cuts, are the minimum quotient cut 
problem and the sparsest cut problem. These ask respectively to minimize the 
cut quotient q(S, S) and the cut density d(S, S) over all cuts [S, S] of a graph G, 
where 

9 ( S) S) = §| and d(S,S)~ l[ ' S '" S ' J 



b(S,S) v ' ' b(S,S)(l-b(S,S))' 
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Both q and d penalise unbalanced cuts, although q does so to a greater extent 
than d. 

Problems such as the minimum quotient cut problem and the sparsest cut 
problem underlie many divide and conquer algorithms [?] , and find applications 
in VLSI layout problems, packet routing in distributed networking, clustering, 
and so on. Unfortunately, for general graphs, finding a minimum quotient cut 
or a sparsest cut is known to be NP-hard 6 ID] . Thus there are two possible 
ways of developing efficient algorithms for these problems: either by considering 
approximation algorithms or by restricting attention to certain graph classes. 
There has been much research done in finding approximation algorithms for 
these problems. Here, we mention only the seminal paper of Leighton and Rao 
[5] giving a polynomial-time 0(log ?i)-approximation algorithm for the minimum 
quotient cut problem, and the significant improvement in the approximation 
factor to 0(y/\ogn) in a paper of Arora, Rao, and Vazirani [2J. On the hard- 
ness side, Ambiihl et al. [1] proved that the sparsest cut problem admits no 
polynomial-time approximation scheme unless NP-hard problems can be solved 
in randomized subexponential time. 

We approach the problems of minimum quotient cut and sparsest cut from 
the perspective of developing exact polynomial-time algorithms for restricted 
graph classes. Such approaches have not received as much attention in recent 
years as the development of approximation algorithms, but we hope this paper 
will take a step towards sparking interest. 

Bonsma [3] gave polynomial-time algorithms for finding sparsest cuts of unit 
circular graphs and cactus graphs. Park and Phillips [13] . building on the work 
of Rao [14] . gave a polynomial-time algorithm for determining the minimum 
quotient cut (as it has been defined above) of planar graphs. Given that many 
planar-graph algorithms have been adapted for generalizations of planar graphs 
- see for example the introduction to [4] and the references therein - surpris- 
ingly little is known about the complexity of computing minimum quotient cuts 
or sparsest cuts for generalizations of planar graphs. Here, we generalize the algo- 
rithm of Park and Phillips to give the first exact polynomial-time algorithm for 
determining minimum quotient cuts and sparsest cuts of bounded-genus graphs. 



1.2 Results 

Before we state our result precisely, we give a generalization of the minimum 
quotient cut and sparsest cut. Notice that the denominators for both the cut 
quotient q and the cut density d are concave and increasing functions of b(S, S) 
on the interval [0, |]. For any concave, increasing function / : [0, |] -> [0, oo 
and a cut [S, S] of a graph G, we define 



db(S,S) 



f{b(S,S))> 
and we define 

d f (G)=mind f G (S, S), 
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where the minimum is taken over all cuts [S, S] of G. Any cut [S, S] that mini- 
mizes d G is referred to as an f -sparsest cut of G. 

Let / : [0, 5] — > [0, 00) be a fixed concave, increasing function that is com- 
putable in polynomial time on the rationals, and let g be a fixed non-negative 
integer. The input for our algorithm is an n-vertex undirected multigraph G of 
genus g. Our algorithm computes an /-sparsest cut of G in time 0(n 2a + 4 f+ 7 ). 



1.3 Overview and Techniques 

In this section we give an informal overview of our methods. Our methods ex- 
tend those of Park and Phillips [13] and combine them with surface homology 
techniques, used for example in [4]. 

A simple averaging argument shows that, given a graph G, there exists a 
sparsest cut [S, S] of G that is minimal, i.e. a cut where the graphs induced by 
G on S and S are both connected. This extends easily to /-sparsest cuts, where 
/ is a concave increasing function. 

There is a standard correspondence between the cuts of a planar graph G 
and the cycles of its dual D(G): the minimal cuts of G correspond precisely 
to the cycles of D(G), and the size of a cut in G is equal to the length of its 
corresponding cycle in D{G). One can similarly construct a dual graph D(G) for 
a graph G embedded on a surface; however the correspondence between cuts of 
G and cycles of D(G) is not quite so simple. Roughly, for a graph G of genus g, a 
cut of G corresponds to a union of at most g + 1 cycles of D(G), but the reverse 
does not hold: a union of at most g + 1 cycles in D(G) does not necessarily 
correspond to a cut in G. Using the surface embedding of G, we construct a 
function from the set of oriented edges of D(G) to I? 9 with the following 
property: summing around the oriented edges of a union of cycles of D(G) 
gives the zero vector if and only if that union of cycles corresponds to a (certain 
generalization of a) cut of G. 

Extending and simplifying an idea from [T3] , we also construct a function w 
from the set of oriented edges of D(G) to Z with the following property: if a 
union of cycles in D(G) corresponds to a cut in G, then summing w around the 
oriented edges of cycles in the union essentially gives the balance of [S, S] . 

Using D{G), 0, and w, we construct a type of covering graph H, again 
extending an idea in [13] . For each fixed value v and k of and w, we can use 
H to find a shortest cycle in D(G) whose (9-value is v and whose w- value is 
fc. Such a shortest cycle of D(G) corresponds to a shortest path in H between 
suitable vertices. 

By repeatedly applying a shortest-path algorithm to H, we obtain, for every 
v and k (in a suitable range), a shortest cycle of D(G) whose 0- value is v and 
whose w- value is k. We construct the set X of every union of at most g + 1 of 
these shortest cycles. The size of X is n°( g \ and we show that at least one 
element of X corresponds to an /-sparsest cut of G. 
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2 Preliminaries 



We begin this section by proving some simple inequalities for concave functions. 
Throughout, rather than working with concave increasing functions / : [0, |] — > 
[0,oo), we work instead with functions / : [0,1] — > [0,oo) that are concave and 
increasing on [0, 5] and have the property that f(x) = f(l — x) for all x G [0, 1]. 
We work with these functions purely for the convenience of having, for any cut 
[S, S] of G, that 

f(b(S,S)) = f(\S\/\V\) = f(\S\/\V\). 

Note that such functions are in fact concave on their entire domain. For the algo- 
rithm, we assume that / can be computed in polynomial time on the rationals. 

Lemma 1. Let f : [0, 1] — > [0, 00) be a concave increasing function on [0, \\ with 
the property that f(x) = /(l — x) for all x G [0, 1] . Suppose X \,...,Xk € [—1, 1] 
and x :=J2 i=1 Xi £ [—1,1]. Then 

k 

/(M)<£/(N). 

i=l 

Proof. We may assume without loss of generality that x G [0, 1] by switching 
the signs of the Xi if necessary. We also know that / is concave on its entire 
domain. Recall that a function / : [0, 1] — > [0, 00) is concave if and only if, for 
every a,b,t <= [0, 1], we have that tf (a) + (1 - t)f{b) < f(ta + (l - t)b). 

We prove the inequality using induction. The case k — 1 is trivial. We prove 
the case k = 2. We have, without loss of generality, the two cases 

(i) xi, X2 G [0, 1] and 

(ii) Xl G [0,1], x 2 G [-1,0]. 

Case (i): Given xi,x 2 G [0,1] with x = x\ + x 2 G [0,1], choose r, s G [0,1] 
such that xi — rx and x 2 = sx; thus r + s = 1. From the concavity of /, we have 

rf(x) + (1 - r)/(0) < f( Xl ) and sf(x) + (1 - s)/(0) < f(x 2 ). 

Adding the two inequalities together, and using the fact that r + s = 1, we obtain 

f(xi) + f(x 2 )>f(x) + f(0)>f(x), 

as required. 

Case (ii): Given x\ G [0, 1], x 2 G [—1,0] with x = x\ + x 2 G [0, 1] apply case 
(i) to the numbers 1 — x\, ~x 2 G [0, 1]. Thus, we have 

f(\xi\) + f(\x 2 \) = f(l - Xl ) + f(-x 2 ) > f(l - Xl - x 2 ) - f(l -x) = fix), 

as required. This proves the case k = 2. 
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The induction step follows easily. Indeed, order the Xi such that 52»=i x i ^ 
[-1,1] for all r = l,...,k. Then 

k fc-1 

/(*) = /(!>) = /(E^+^) 

i=l i=l 
fe-1 

< /(|E^|) +/(kfc|) (by case k = 2) 

i=l 

fc 

< 2_, f(\ x i\)- (by induction hypothesis) 

z=l 

□ 

Next we prove that for any graph G, d G can be minimized by a minimal 
cut, that is, a cut [5, «S] for which both G[S] and G[S] are connected (here G[.A] 
denotes the graph induced by G on A C K). This is a well-known fact for the 
cut quotient q and the cut density d, and generalises easily to d G . Before we can 
do this, we need a trivial averaging argument; we state it formally so that we 
can refer to it later. 

Proposition 1. For i = 1, . . . , k, let ai, pi, and qi be non-negative real numbers 
with at least one qi non-zero. Then there exists some %' G {l,...,fc} such that 
qv ^ and 

VV_ < EjU a iPi 
Qi' ~ Ej=l 

Proof. Choose i' to minimize Pi'/qe (if qi = 0, we take Pi/qt to be oo). Thus 
we have that Pi>/qi> < pi/qi for all i = l,...,k. Rearranging, we have that 
Pi'(a,iqi) < (aipi)qi> for all i = 1, . . . , fc. Summing both sides over i and rearrang- 
ing gives the desired inequality. □ 

Proposition 2. Let G = (V, E) be a connected graph and let f : [0, 1] — > [0, oo) 
be as in Lemma[l[ Then there exists a cut [S, S] of G such that dg(5, S) = d? (G) 
and for which G[S] and G[S] are both connected. 

Proof. Amongst all cuts of G minimizing d^, let [A, B] be one of minimum 
size. Suppose that G[A] has components A\, . . . , A r and G[B] has components 
B\ , . . . , B s . We claim that there is some i' for which d f G {A v ,Ai.) < d f G {A,B) 
(and by symmetry, there is some j' for which d G (Bjt, Bji) < d G (A, B)). Since G 
is connected and [A, B] is of minimum size, the claim implies that both G[A] and 
G[B] have only a single component, proving the proposition. It remains to prove 
the claim. Let \V\ — n, ai — \Ai\/n, and bi = \Bi\jn. Let a — \A\/n — J2l=i a i 
and b = \B\/n = J2t=i Let 

e« = \{ab € E ; a e A h b e BA\ and d{, = 

}{ai) 
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Wc have that 



d f (G) = dUA,B) = l [A ' B] \ = 

{ ' G{ ' f(b(A,B)) /(EI=i«i) 



ELi /K) E^=i 4 



/(E[=i a O 

> 



ELi /(«<) EJ=i 4- 



EU/K) ' 

where the inequality follows from Lemma [T] Applying Proposition [1] we find 
that there exists some i' such that 



/(Oi') /(Oi') 



□ 



It turns out that the description and the proof of correctness of our algorithm 
is most conveniently and naturally expressed in the language of surface homology. 
Through the remainder of this section, we introduce the necessary concepts 
keeping our treatment as simple and self-contained as possible. One can find 
more comprehensive treatments in e.g. [7|8j . 

Although we are only concerned with undirected graphs when determining 
quotient cuts, sparsest cuts, and other vulnerability measures, we shall have 
cause to orient edges of our graph through the course of our proofs and algo- 
rithms. Each edge e = ab of a graph G = (V, E) has two orientations, namely 
(a,e,b) and (b, e, a). The two orientations are denoted it and V, although we 
cannot say which is which in general. Given a set E of edges, we write E for 
the set of their orientations, two for each edge. The edge space of G, denoted 
£(G), is the free abelian group on E modulo the relation that ^e" = —if. For 
each p £ £(G), we can express p uniquely as 

where A-^ G Z for all ~£ € E and min(A-^-, A%) = 0. Then, we define 

Thus \p\ in a sense counts the number of edges in p. 

The cut space T{G) of G = (V,E), which is a subgroup of £ (G), is defined 
as follows. For a cut [S, S] of G, define 

3 \ - 

[S,S] = 2^ (a,e,b). 

ab=e£E 
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Then T(G) is the subgroup of 8(G) generated by {[S,5]:5C V}. 

The next lemma shows how, by suitably assigning weights to oriented edges 
of a graph, we can determine the balance of a cut simply by summing the weights 
of the edges in the (oriented) cut. This is a generalisation of a result for planar 
graphs that was presented (but not proved) in [13]. 

Lemma 2. Let G = (V, E) be a connected graph and let v G V be some fixed 
vertex. There exists a function w : £^ — > Z with the following properties. 

(i) For every ~e G ~e] , we have w(^e) = —w(~e). 
(ii) For every ~e( E ~e" , we have \w(~e*)\ < \V\. 
(Hi) For every S satisfying v € S C V , we have that 



Furthermore, a function satisfying the above properties can be constructed in 
0(n 2 ) time. 

Remark 1. The function w described in Lemma [2] can be extended to a homo- 
morphism w : £{G) — > Z because of property (i). 

Proof. Let T = (V,Et) be a spanning tree of G and let v be a root of T. If the 
lemma holds for T, then it clearly holds for G by setting w(a, e,b) =0 for all 
ab e E\E T . 

Let ab £ Et, and without loss of generality, assume that b is a descendant 
of a (that is b is further from v than a). Deleting ab disconnects T into two 
components S a b and <§ a f,, where S a b is the component not containing v (and 
hence not containing a). Set w(a, e, b) = —w(b,e,a) = \S a i>\. We do this for 
every ab G Et- 

Clearly w satisfies properties (i) and (ii), and furthermore, it is not hard to 
see that w can be constructed in 0(n 2 ) time. It remains only to prove property 
(iii). We prove the claim by induction. Let v\, . . . , v r be the vertices of T adjacent 
to v. Let Ti — (Vi,Ei) be the subtree of T formed from vi and its descendants 
for i = 1, . . . ,r. Given v € S C V and S = V\S, let Si = SnVi and $ = Sn Vi 
Also, let 



w([S,S}):= ]T u;(o,e,6) = |5|. 



(a,e,b)£[S,S} 



M = {i : tn e S}. 



Observe that 



r 




i=l 
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Now we have 



w([S, S]) = y]w([Si, Sj]) + 2J w(v,Wi,Vi) 

i=l i^Af 

= X! 1^*1 ~~ E 1^1 + E w ( v , vv i, v i) (induction) 

ieM ig_M i<£M 



£N + E(l y *l-N) 

r 

£i&i = isi. 



□ 



Throughout, w will be a homomorphism from £(G) to Z satisfying the prop- 
erties of Lemma [5J 

The domain for the function d G is the set of all cuts of G. It turns out that 
it is necessary to extend d G in a natural way to a function on T(G). Minimizing 
d G on T(G) will be equivalent to minimizing it on the set of cuts of G. Although 
T(G) is an infinite set, we shall eventually look to minimize d G over a suitable 
finite subset of T(G) . 

For each <j> e T(G), if w(<fr) ^ 0, we define 



where n = \V\; if f(\w(<f>)\) = 0, we define d f G (<f>) = 00. Note that if <f> = [S,S], 
then \(j>\ — |[5, 5]| and |io(^)| is either |5| or |5|. Hence the function d G defined 
above extends the definition of d G given in the introduction. 

Our next lemma shows that minimizing d G over T(G) is equivalent to min- 
imizing d G over simple cuts [5, S] of G. First some observations. For the cut 
space T(G) of G = (V,E), we note that if 5 is the union of disjoint subsets 5i 
and S2 of V, then 

\0\ = [s u Si] + [s 2 ,si 



Note also that [5, 5] = — [5, S]. Thus every element of T(G) can be written as a 
positive integer linear combination of single vertex cuts [{v}, {v}]. 

Lemma 3. Let G = (V, E) be a graph. For every 4> G T{G), there exists 5 C V 
such that 

4(5,5) < 4 (</>), 

so in particular 

min 4(0) = d f (G). 
4>eT(G) GV ' K ' 
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The proof uses the same averaging argument used in the proof of Proposition^ 

Proof. The lemma clearly holds if d G {(j)) — oo, so assume <f(j(^>) < oo. Choose 
an integer fc such that 

where < A„ < fc for all v E V. For each i — 1, . . . , fc, let 



Uw- 



Thus we have 5i 3 52 3 ■ ■ ■ 3 ft, and 

fc 



i=l 

Since the sets are nested, there is no cancelling of edges; hence 

fc 



|0|=£|[^I]|. 



i=l 

Using the fact that w is a homomorphism together with the triangle inequality, 
we have 

fc fc 



w 



,(cj>)\ = \j2H[Si,Si])\<J2\ w ti s i>^ 

i=l 

Now we have 



d fu) = M > 



> ^tilt-g^jl . Lemma □ 

Etr/d^^l) 



^ — I [ y 4 ] _l — gome ji Proposition [T] 

" /( r ([^]) |) 

□ 

We specify a roott u; of a graph G = (V, E) by giving an alternating sequence 
of vertices and edges w = (xx, e\,X2, e2, . . . , %k-i, efc-i, £fc), where (x*, e^, x,+i) G 
1^ for all i. We write |u>| for the number of edges traversed in w (which in this 
case is k — 1). H xx = Xk then w is called a closed walk. If all edges of a closed 
walk w are distinct, then w is called a called a circuit of G. If all the vertices of 
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a closed walk w are distinct (except x\ = Xk), then w is called a cycle. A walk 
in which all vertices are distinct is called a path. We write ~v$ for the element of 
£{G) given by 

fc-i 

wt = ^(x 4 ,e 4 ,x l+ i); 

i=l 

we refer to as an oriented walk, (circuit, etc). Note that for a walk ty, we have 
| to | > | with equality when w is a circuit. 

The cycle space C(G) is the subgroup of £(G) (redundantly) generated by 
the oriented cycles c of G. Note that C(G) contains all oriented closed walks of 
G. 

We now turn our attention to graphs embedded on closed orientable surfaces. 
Formally, a surface is a compact connected topological space in which every 
point of the surface has an open neighbourhood homeomorphic to R 2 or the 
closed halfplane {(x,y) € R 2 : y > 0}. The set of points having halfplane open 
neighbourhoods is called the boundary of the surface. Every component of the 
boundary is homeomorphic to the circle S . A closed surface is one without 
boundary. A surface is called orientable if it does not contain a subset (with the 
subset topology) homeomorphic to the mdbius band. 

The genus of a connected, orientable surface is the maximum number of 
cuttings along non-intersecting closed simple curves that can be made without 
disconnecting the surface. It is well known from the classification of surfaces that 
every closed orientable surface of genus g is homeomorphic to a sphere with g 
handles. Thus we can think of every such surface as embedded in R 3 . Informally, 
a graph G can be embedded on a surface E if G can be drawn on E in such 
a way that no edge crosses a vertex or another edge, except possibly at its end 
points. For example, all planar graphs can be embedded on the sphere. 

More formally, we have the following definitions. A (topological) path in E 
is a continuous function 7 : [0, 1] — > E, where 7(0) and 7(1) are called the end 
points of the path. If 7(0) = 7(1) then 7 is called a (topological) cycle. Paths and 
cycles are referred to collectively as curves. A curve is simple if it is injective, 
except in the case of a cycle where we permit 7(0) = 7(1). We often do not 
distinguish between curves and their images in E. 

Throughout, we shall consider embeddings of multigraphs with loops on ori- 
entable surfaces. (Although G has no loops, its dual graph may have loops as 
we shall discuss later.) An embedding of a graph G in E maps vertices v to dis- 
tinct points ip(v) of E and maps edges (resp. loops) e = ab to simple topological 
paths (resp. cycles) -f e of E: two such paths may intersect (if at all) only at 
their end points, and the end points of -f e are precisely tp(a) and ip(b). Oriented 
edges if — (a, e, b) are embedded by simple paths 7-^, where we insist that 
7-^(0) = ip(a) and 7-^(1) = "4>{b). We often abuse terminology and notation by 
identifying vertices, edges, and walks of G with their images in the embedding 
of G on E. 

From an algorithmic point of view, one can use rotation systems to input or 
output embeddings of graphs on surfaces. We do not define rotation systems here 
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because we shall only use graph embeddings indirectly when applying existing 
algorithms to our problem; instead we refer the reader to [T2] . 

Throughout, we shall only consider cellular embeddings. An embedding of 
G on £ is called cellular if removing the image of G from £ leaves a set of 
topological disks called the faces of G. The genus of a graph G is defined to 
be the smallest integer g such that G can be embedded on a closed orientable 
surface of genus g. If G is a graph of genus g, then every embedding of G on a 
closed orientable surface of genus g is cellular (Proposition 3.4.1 [12 ), and for 
fixed g, such an embedding can be found in linear time [llj . Euler's Theorem 
gives the following relationship between the number of vertices n, the number 
of edges m, the number of faces £, and the number of boundary components b 
in a cellular embedding of a graph G on a surface £ of genus g: 

n-m + £ + b = 2-2g. 

We now define the boundary space of a graph G embedded on £ . Each ori- 
ented edge of G separates two (possibly equal) faces of G denoted left(e) 
and right (it). (The notion of left and right with respect to an oriented edge 
is well defined for orientable surfaces.) For a face F, the oriented edges e for 
which left(e) = F taken in order (with appropriate intervenng vertices) form 
a closed walk / around F, which we call the facial walk of F. Let Fx, . . . ,Fg be 
the faces of the embedding and let fi be the facial walk of Fi . The oriented facial 
walk fi of Fi is given by 

t= £ 

~t: left(~et)=Fi 

Notice that /{ may contain two oppositely oriented edges, but such edges cancel 
in fi , leaving the sum of oriented edges that form the boundary of Fi . The 
boundary space B(G, £) is defined to be the group generated by fx, ■ ■ ■ , ft and 
is easily seen to be a subgroup of C(G). Note that the boundary space of G, in 
contrast to the cycle space and cut space, depends on the embedding of G. 

The geometric dual of a graph G cellularly embedded on £ is denoted by 
D(G) = (V',E') and is the graph (with cellular embedding on £) constructed 
from G as follows. For each face F of G, a vertex D(F) is placed inside F: 
these are the vertices of D(G). Each edge e of G has a corresponding edge D(e) 
in D(G): D(e) = D(Fx)D(F2), where e is the edge separating the (possibly 
indistinct) faces Fx and F2 (and D(e) crosses e and no other edge of G in the 
embedding of D(G)). Note that G (with its embedding) is a dual of D(G) (with 
its embedding). Therefore, each vertex v of G corresponds to a face D(v) of 
D(G). Thus D maps vertices, edges, and faces of G bijectively to faces, edges, 
and vertices of D(G) respectively. We extend D to map oriented edges of G to 
oriented edges of D(G) as follows. Given an oriented edge e , we set D(e) = 
(D(left(~ef)),D(e),D(right(~ef))). This, however, reverses the sense of left and 
right so that D(D(lt)) = -it. 

Remark 2. Although G is a loopless graph, D(G) may not be. Nonetheless, all 
notions introduced so far carry through naturally for loops of embedded graphs. 
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In particular, a loop e of an embedded graph G has two orientations it and V, 
although we cannot specify which is which by the order of the end vertices. If 
D(e) is a loop for some edge e of G, then D(e) should cross e from left to 
right. 

Since D bijectively maps edges of G to edges of D(G), we see that D can 
be extended to an isomorphism D : £{G) — > £(D{G)). We have the following 
well-known correspondence. 

Theorem 1. The restriction of D to T{G) gives an isomorphism T{G) — > 
B(D(G)). 

Proof. Note that for a vertex v of G, and f v the facial walk around D(v), we 
have 

0([M, {«}]) = -£■ 

This defines a bijective correspondence. □ 

Rather than working with T(G), we can work instead with B(D{G)) using 
the isomorphism I?. Since all boundaries are sums of oriented cycles, we can 
use shortest-path algorithms to find shortest boundaries in D{G), which if done 
suitably, can give us /-sparsest cuts in G. 

For a e B(D(G)), define w(a) := w(L»- 1 (cr)). Notice that for all a e 
B(D(G)), we have |-D"V)I = |ct|. Thus, defining 

*^ (ff) = ?(i&)' (2) 

we have that 

min dL.^Aa) = min dt,(d>) = d f (G). 

a<LB(D(G)) D< - G > y ' 0GT(G) GK ' ' 

We now set about trying to minimize d? D , G y Our next lemma says that when 
minimizing $ D , G y we can restrict attention to elements of B{D(G)) that are the 
sum of at most g + 1 oriented circuits (where g is the genus of G) . 

First a proposition. 

Proposition 3. Let G(V, E) be a graph cellularly embedded on E. For every cut 
[S 1 , S] of G there exist vertex disjoint circuits w\, . . . , w r such that 



a = D([S,S\) = w[ + ---+- 
and 

\a\ = \wt\ H h \wt\. 



Proof. Proof. We know a = D([S,S}) is an element of B(D{G)) C C(D(G)) in 
which no edge occurs more than once. Thus we can write 

s 

1=1 itet 
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where each cf is an oriented cycle and each A-^ G {0, 1}. The number of oriented 
edges of a entering and exiting any given vertex must be equal since this is the 
case for any oriented cycle and remains the case after cancellation of oriented 
edges when summing oriented cycles. Thus we can walk around D(G) using the 
oriented edges of {if : A-g> = 1} to form disjoint closed circuits w\, . . . , w r , where 
a = w\ + ■ ■ ■ + w^. Furthermore, there is no cancellation of edges when we add 
these disjoint oriented circuits together; hence 

M = |^| + --- + |^|. 

□ 

Recall that for each vertex v of G, D{y) is a face (a topological disk in E) of 
D(G) and f v is the facial walk of D(v). The closed face D*(v) is a closed disk 
which has the facial walk /„ embedded along its boundary. Every G ~~E occurs 
on the boundary of some D*(v); thus each edge of G is either embedded on two 
distinct closed faces or is embedded twice on the same closed face. Note that E 
can be constructed by gluing these closed faces together along common edges of 
G (respecting the orientation). We can now prove our lemma. 

Lemma 4. Suppose G — (V, E) is a graph cellularly embedded on a surface E 
of genus g. There exists a G B(D(G)) that minimizes d, D r G <* such that 

a = wi H h w r , 

where w\, . . . , w$ are disjoint oriented circuits of D(G) and r < g + 1. Further- 
more m > \a\ = \w\\ + ■ ■ ■ + \wl\, where m = \E\. 

Proof. By Proposition [21 there exists an /-sparsest cut [S, S] of G for which 
G[S] and G[S] are both connected. Thus [S, S] G T(G) minimizes d G (over all 
elements of T(G) by Lemma H, and so a — D([S, S]) minimizes d D ^ (over 
all elements of B(D(G))). By Proposition [31 we know that we can write a = 

D([S, S]) = wt + ■ ■ ■ + w^, where m > \a\ = |toi| + ■ ■ • + \wt\. It remains to show 
that r < g + 1. 

Let (Hs, Ss) be obtained from (G, S) as follows. Take the set of closed faces 
{D*(v) : v G S} and for each edge v\v 2 G G[S], glue D*(v\) and D*{v 2 ) together 
along their common edge D{v\V2), respecting orientation, to form E$- Since G[S] 
is connected, Us is a surface with boundary. The facial cycles of the glued faces 
now give a graph Hs (which is a subgraph of D{G)) embedded on Es- Observe 
that each edge of D(G) occurs at most once as an edge of Hs and the edges of 
wi, . . . ,w r form the boundary of Es- Disjoint edges of D(G), when they occur 
in Hs, remain disjoint; hence, since w\, . . . ,w r are disjoint, Es must have at 
least r boundary components. We form {Hs, Eg) analogously and by symmetry 
it has the same properties described above. Note that E can be formed from 
Es and E§ by gluing them together suitably along their boundaries. Note also 
that the edges of Hs are precisely the duals of the edges of G incident with S: 
similarly for H§- We apply Euler's theorem to the two embeddings. 
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Let <?, gs and gg he the genii of E, E$, and Eg respectively. Let n, m, £, and 
6 = denote the number of vertices, edges, faces, and boundary components 
for the embedding of G on £. Let ns, nis, is, and b$ denote the numbers of 
vertices, edges, faces, and boundary components respectively for the embedding 
of Hs on Us, and analogously for S. Finally, let rib and mb denote the numbers 
of vertices and edges that occur on the boundaries of Eg and Eg. We have 



Using Euler's formula and the first three equalities, we have 

2-2g = n-m + £ + b=(2- 2g s ) + (2 - 2g s ) + m b - n b - b s - b§ 
Rearranging, and using the fourth and fifth statements above, we have 



We can easily use shortest-path algorithms to find shortest cycles in a graph. 
However in this situation, we are required to find a shortest boundary (loosely 
speaking). We require a simple way of testing whether a cycle is a boundary. This 
is accomplished using the ideas of homology. We only require the following fact 
about the homology of graphs on surfaces. If G is a graph cellularly embedded 
on an orientable surface E of genus g, then the quotient group C(G)/B(G) is 
isomorphic to I? 9 ; this follows easily from standard results on cellular homology. 

For an n- vertex graph G cellularly embedded on a closed orientable surface E 
of genus g, a system of loops is a set of cycles of G through a common vertex such 
that cutting E along these cycles gives a topological disk. By Euler's formula, 
every system of loops must consist of 2g cycles. Erickson and Whittlesey [5] give 
a greedy algorithm, which, given a cellular embedding of G on E, finds a system 
of loops c(l), . . . , c(2g) in 0(n log n + gn) time. 

We define a homomorphism @, : C(D(G)) — > Z, where, for every a £ 
C(D(G)), the integer 0i(a) measures the net number of times a crosses c(i) 
(here sign indicates the direction of crossings). Let us define ©i formally. 

For each oriented edge if of G and each oriented edge er of D(G), define 



n = rii + n 2 - n bl 
m = mi + TO2 — mb, 

e = e! + e 2 , 

rib = m,b (since boundaries consist of disjoint cycles) 



bs, bg>r 



S + 1 = 9s + 9s + 2 ~ mfc ) + + b s) J ^ r - 



□ 




1 



if D(Y) 
1 i£D(t) 








otherwise. 
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For every <fr — J2i *= anc ^ everv c = J2j Mj^f e £(D(G)), we define 



which counts the directed number of times </> and 7r cross each other. It is easy 
to check that the above is well defined. We write 0i as a shorthand for 

c(») 

Finally, define : £{D{G)) -> Z 2 s by setting 

0(a) = (0 1 (<7),...,0 2s (a)) 

for every cr e £{D{G)). 

We have the following proposition which effectively says that any cycle of 
D(G) that intersects each c(i) a net number of zero times is a boundary. It is 
very much what we expect from the properties of homology, but we give the 
details for completeness. 

Proposition 4. The function defined above is a surjective homomorphism 
from C(D(G)) to Z 2 » whose kernel is B(D(G)). 

Proof. The fact that is well defined and a homomorphism is easy to check. 

To see that is surjective, we must find for each j = l,...,2g, a cycle 
c*(j) € C(D(G)) such that 0(c*(j)) = ±Uj, where u, is the vector that has 1 
in the jth component and O's elsewhere. Consider the embedding of G on the 
(closed) topological disk T formed by ungluing S along the cycles c(l), . . . , c(2g); 
thus each unglued edge is embedded twice along the boundary of T. Note that 
each c(j) has at least one edge ej that does not belong to any of the other 
c(i)'s: indeed, if all the edges of c(j) belonged to other cycles, then c(j) would 
be redundant and we could cut the surface into a topological disk with fewer 
than 2g cycles. 

Pick points x and x on the boundary of T that lie in the interior of the 
two embeddings of ej. Let 7 be a topological path in T from x to x that is not 
incident with any other points along the boundary of T and not incident with 
any vertices of G. It is clear that such a path exists. Then 7 corresponds to a 
topological cycle in S that crosses the cycle c(j) exactly once and does not cross 
any of the other c(i)'s. By listing the alternating sequence of faces and edges of 
G that 7 crosses, we obtain a cycle c*(j) in D(G) and hence an oriented cycle 

c*(j) € C{D(G)) that crosses c(j) once but does not cross any of the other c(i)'s. 
Thus, we have 



0(c*(j))=±u j , 

showing that is surjective. 

Finally we show that ker(6>) = B(D(G)). For any oriented cycle and any 
oriented facial walk ~f~ of D(G), we have 
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since ~C* enters and exits the face Fi an equal number of times so that all contri- 
butions cancel. Thus, we have that 0(j) = for all oriented facial walks 

7 of 

D(G). Therefore 0(<j>) = for all G B(D(G)) showing that B(D(G)) C ker(<9). 
Now we have 

Z 2 3 s C(D(G)) ^ C(D(G))/B(D(G)) ^ J?* 



ker(</>) ker(^)/B(D(G)) ker(0)/B(£>(G)) ' 

where the first isomorphism is from the First Isomorphism Theorem, the second 
is from the Third Isomorphism Theorem, and the third is given by the fact 
about homology groups given earlier. We deduce that ker (<fi) / B(D(G)) must be 
the trivial group (by standard properties of finitely generated abelian groups) 
and that therefore we must have ker(0) = B(D(G)). □ 

For convenience, we combine some of the results we have so far to give the 
following proposition. 

Proposition 5. There exists an element a £ B(D(G)) that minimizes d? D , G -, 
and satisfies the following properties. 

(a) We can write a as W\ + • • • + wt where the w\ are oriented circuits in D(G) 
and r < g + 1 . 

(b) We have \o~\ = \wi\ + • • ■ + \w r \ with \wi\ < m for all i. 

(c) 0(a) = 0{w\) + ■ ■ ■ + 0{vt) = 0. 

Proof. Statements (a) and (b) follow from Lemma 2J and statement (c) follows 
from Proposition |4j □ 

Next we define a covering graph of D(G) in which certain shortest paths will 
correspond to elements of C(D(G)) whose sum will minimize d D , G y 

We construct an infinite multigraph H = (Vh,Eh) from D(G) as follows. 
We set Vjj = V x Z x I? 9 , where V is the vertex set of D(G). For convenience, 
we describe the oriented edges of H before describing its edges. For each oriented 
edge e = (u%,e, u^) of D(G) and each (fc,v) G Z x Z 2s , we have an oriented 
edge of H, denoted (it, k, v)*, from 

(ui,fc,v) to (u2,k + w( e),v + @( e)). 

The same edge oriented oppositely is given by ( i e',k + w('^),y + 0(~e ¥ ))*. We 
write ( e , k, v) for the edge of H corresponding to the oriented edge (it, k, v)*; 
thus every edge of H has two labels. 

Walks of D(G) naturally correspond to walks of H as follows. Let w — 
(uq, ei,u\,..., et-i, ut) be a walk of D(G). Let it;, = (uq, e\, tti, . . . , Ui) be 
the same walk up to the ith vertex, and let et = (wi-i, ej, Ui). Let H(w) be the 
walk in H given by H(w) = (uq , , uf , . . . , e$r_i, ), where Uq = (u, 0,0) 
and 

uf = (ui,w(Wi),Q(wi)) and ef = (et,w(wi-i),0(wi-l)). 
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It is easy to check that w i— > H(w) is a bijective correspondence between walks 
of D(G) and walks of if that start at (tt, 0,0) for some u E V . Furthermore, w 
is a walk of D(G) from u to v! and satisfies u>(wt) = k and @(io) = v if and 
only if H(w) is a walk of if from (u, 0, 0) to (r/, k, v). 

Defining Vj*- = V x {— mn, . . . , mn} x {— m, . . . , m} 29 C Vh , let if* be the 
finite graph induced by if on V^. For (w, fc, v) £ V^, define p(u, k, v) to be a 
shortest path in if* from (u, 0, 0) to (u, k, v) (if it exists); this path corresponds 
to a closed walk in D(G). For fixed k and v, let p(k, v) be the path of minimum 
length in {p(u, k, v) : u € £/} (if it exists). Let w(k,v) be the (closed) walk of 
D(G) corresponding to the path p(k, v) in if*. We have the following lemma. 

Lemma 5. Suppose w is a circuit of D(G) such that w(w) = k and &(ln) = v 
and 1 10 1 < to. TTien 

|t£| > |w(fc,vj|. 

Proof. Suppose the circuit w starts at u € U. We shall show that H(w) is a walk 
on if* (from (u, 0,0) to (u,fc,v)). Then using the fact that if is a circuit and 
consequently |to| = \w\, we have 

|t*| = |w| = \H(w)\ > \p(u,k,v)\ > \p(k,v)\ = \w(k,v)\ > \w(k,v)\, 

proving the lemma. 

To see that H(w) is a walk on if*, we note first that \~v$\ — \w\ < m. Let Wi 
be the same walk as w up to the ith vertex. Then by property (ii) from LemmaEJ 
we have that |m)(w^)| < \wi\n < mn. Also, since ~w < m, it is not hard to see 
that G(wi) G {—m, . . . ,m} 29 . This shows that the walk if(io) never leaves V^, 
as required. 

□ 



Let W be the set of all the w(k, v). Let 

X = {w|h ^wt:WieWVi and r<g+l}, 

and let Y = {a € X : 0(a) — 0}. We have the following corollary. 
Corollary 1. There exists an element ofY that minimizes dp( G y 

Proof. Clearly every element of Y is in B(D(G)). Let a £ B(D(G)) minimize 
d(G) anc ^ sa tisfy the three conditions of Proposition [SJ in particular 

a = Wi H h Wr 

for r < g + 1, where each < to, and where O(o~) = 0. For each i, let 
6*(wi) = v; and let w(wi) = ki. Then by Lemma[5]we have that 

\w(ki,Vi)\ < \v%\. 
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Let a* £ X he given by 

a* = w(ki,vi) + • • • + w(k r ,v r ). 

We have that 

\<r\ = \wx\ + ■■• + \wt\ > |w(fci,vi)| + --- + |u;(fc r ,v r )| > |ct*|, 
that 0(a*) = 0(a) = 0, and that u>(<7*) = w(a). Thus a* G Y and 

showing that tr* € Y" minimizes d G . □ 

We now have all the ingredients to present our algorithm and to prove its 
correctness. 

3 The Algorithm 

In this section, we present the basic steps of our algorithm and compute its 
running time. In order to keep the presentation simple, we do not optimize the 
running time. Our algorithm runs in time 0(n 2g + 4 s+ 7 ) ) and it seems unlikely 
that our methods can give an n°^ 9 )-time algorithm without significant modifi- 
cation. 

Let / : [0, |] — > [0, oo) be a fixed concave, increasing function that is com- 
putable in polynomial time on the rationals, and let g be a fixed non-negative 
integer. The input for our algorithm is an n-vertex undirected multigraph G of 
genus g. Since n = \V\ then m = \E\ = 0(n) from Euler's formula. Using the 
result of Mohar [11] mentioned earlier, we can find an embedding of G on a 
surface S of genus g in 0(n) time. From the embedding we can construct (in 
0(n) time) the dual graph D(G) — (V, E') together with the function D which 

maps each oriented edge ~~£ G E of G to its dual D(t) G E' in D(G). We have 
\E'\ = \E\ = 0(n), and from Euler's formula, we have \V'\ = 0(n). By the result 
of Erickson and Whittlesey [5] mentioned earlier, we can find a system of loops 

c(l), . . . , c(2g) of G in time 0(n), At this point the algorithm no longer requires 
the embedding of G. ^ 

Next we construct and store the (restricted) functions w : E' — > Z and 

& : E' —> J? 9 . The computation and storage of these functions imposes an 
insignificant time cost in the final analysis, so any crude bound on the running 
time is sufficient. Recall that a function w : E — > Z satisfying the properties 
of Lemma [5] can be constructed in 0(n 2 ) time. Then w — w o D^ 1 can also 
be constructed in 0(n 2 ) time. In order to construct 0, observe that we can 

construct Q-g : E' — >• Z in 0(n) time. Now Oi can be computed in 0(n 2 ) time 
since if c(i) = e% + • • • + et (where k < n, then 
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Thus = (0i, . . . , 02g) can be computed in 2gO(n 2 ) = 0(n 2 ) time. 

From the (restricted) functions w and 0, we can construct the graph H* 
directly from its definition (given in the previous section). Observe that H* has 
0(n)-(2mn+l)-(2m+l) 2g = 0(n 29+3 ) vertices and 0(n 29+3 ) edges. Thus it takes 
0(n 2g+3 ) time to construct H* . For each (u, k, v) 6 Vh* , we compute and store 
the shortest path p(u, k, v) from (u, 0, 0) to (u, k, v) in H* . Finding each shortest 
path requires 0(\ Vh* | log(|Vjj* |)) = 0(n 29+A ) time using Dijkstra's algorithm, 
and so, computing all the p(u,k,v) requires | Vh* \0(n 29+A ) = 0(n i9+J ) time. 
For each fixed k E {—mn, . . . , run} and v E {— n, . . . , n} 29 , we compute and 
store p(k,v), the path of minimum length amongst the p(u,k,v), and we use 

p(k, v) to compute and store w(k, v) and w(k, v) (recall that w(k, v) is the closed 
walk in D(G) corresponding to p(k, v) as described in the previous section). The 
time cost so far is 0(n A9+1 ). 

Recall the sets W, X, and Y from the previous section. Having stored the 
set W of all walks w(k, v), we compute and store the set 

X = {w|h h w$ : Wi e W Vi and r<g + l}. 

Adding elements of W together requires 0(n) time; hence computing and storing 
X requires 0(n\X\) = 0(n\W\ 9+1 ) = 0(n ■ n 2< - 9+1 ^ 2 ) time. We compute 0(a) 
for every a E X and store the set Y = {a E X : 0(a) = 0}, which takes 
0(n|X|) time. Finally we find an element a* of Y that minimizes d D ^ G y which 
takes 0(n|Y|) time, and from Corollary Q] we know a* minimizes d D , G -> over all 
elements of B(D(G). Thus (G) is given by ^(G)( <7 *)' an< ^ tne t° ta l time taken 
to find a* is dominated by max(0(n 2 (f +1 ) 2+1 ), 0(n 49 + 7 )) = 0(n 2 f 2+4 9+ 7 ). 

In order to find an /-sparest cut of G, we compute <fi* = D^ 1 (a*) and 
decompose it as in the proof of Lemma [3] to give 



i=l 

where k = 0(mn) from our bound on H*. Now one of the cuts [Si, Si] is an 
/-sparsest cut by the proof of Lemma [31 which can be found in 0(n 2 ) time once 
a* has been found. 

4 Open Problems 

An obvious question that arises from this work is whether the running time of 
our algorithm can be improved. Specifically, it would be interesting to know if 
the problem of finding the edge expansion of a graph is fixed parameter tractable 
with respect to genus. 

Our algorithm crucially relies on our graph being embedded on an orientable 
surface. In particular, we use the fact that a graph embedded on an orientable 
surface has a directed dual; this is not the case for graphs embedded on non- 
orientable surfaces. It would be interesting to develop methods for finding edge 
expansion of graphs embedded on non-orientable surfaces. 
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